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Abstract. Quantization of electrodynamics in curved space-time in the Lorenz gauge 
and with arbitrary gauge parameter makes it necessary to study Green functions of non- 
minimal operators with variable coefficients. Starting from the integral representation of 
photon Green functions, we link them to the evaluation of integrals involving T-functions. 
Eventually, the full asymptotic expansion of the Feynman photon Green function at small 
values of the world function, as well as its explicit dependence on the gauge parameter, 
are obtained without adding by hand a mass term to the Faddeev-Popov Lagrangian. 
Coincidence limits of second covariant derivatives of the associated Hadamard function are 
also evaluated, as a first step towards the energy-momentum tensor in the non-minimal 
case. 
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1. Introduction 



The investigation of Green functions has led, over many decades, to several developments 
in quantum field theory [1-8]. In particular, if flat Minkowski space-time is studied, we 
can say in modern language that the bare photon Green function (or bare propagator) is 
obtained as follows: 

(i) Work out the invertible gauge- field operator P^ v acting on the potential A u (x), when 
the action functional consisting of Maxwell term plus gauge-fixing and ghost terms is 
considered in the path integral. 

(ii) Obtain the symbol E(P MI/ ) = T^^ik) of P^ u by working in momentum space. The 
symbol is a 4 x 4 matrix of functions on the cotangent bundle of R 4 . To each <9 M term 
in the original operator there corresponds i/c M in the symbol, while multiplicative parts in 
P^ u remain unaffected. 

(iii) Invert the symbol E M „(fc) to find the matrix T, v \k) such that E^E^ = E^E^ = 5*. 

(iv) Obtain the photon Green function 



where the choice of contour 7 reflects the choice of boundary conditions. 

In curved space-time, however, the gauge-field operator is no longer a constant- 
coefficient partial differential operator, and hence only a local momentum-space represen- 
tation can be achieved, after using Riemann normal coordinates [7] . A valuable alternative 
tool is instead provided by the space-time covariant approach of DeWitt [2] to studying 
field theories possessing infinite-dimensional invariance groups. With modern language, 
our starting point is therefore an action functional S consisting of the Maxwell term, 
plus gauge-fixing term in the Lorenz gauge* plus ghost-field contribution, i.e. (hereafter 

* this is due to L. Lorenz [9], not H. Lorentz 




(1.1) 
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g = -det g^) 



S 



d^x^fg 
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4^ ^ v v pP 2a ^ 



(1.2) 



where = V M A^ — V^A M = d^A^ — d^A^, V M is the covariant derivative with respect to 
the Levi-Civita connection, a is a gauge parameter, x an d V> are independent ghost fields 
obeying Fermi statistics [10]. More generally, the gauge-fixing term [more appropriately 
called gauge averaging, from the point of view of the path integral] for Maxwell theory 

might be written in the form — ' with $ any functional on the space of (gauge) con- 
nection 1-forms A^dx^ such that the resulting gauge-field operator P^ v on A v is invertible, 
but we here choose to work in the Lorenz gauge, following Endo. The gauge parameter 
a also occurs in the ghost action since, following Nielsen and van Nieuwenhuizen [11], to 
the ct-dependent gauge-breaking term there corresponds a non-trivial ct-dependence of the 
ghost effective action. The action (1.2) is invariant under the BRST transformations for 
Abelian theory [8,12] and can be eventually cast in the form 



S 



= J d 4 x^/g 



1 



X 



-A^(a)A„ + ^=P ^ 



a 



(1.3) 



where the second-order differential operators P' ll/ (a) and Pq read as 



(1.4) 



□ ee-^V^V, 



(1.5) 



and we assume that they have no zero-modes. Following the Fock-Schwinger-DeWitt 
method [2,13,14], we now consider two abstract Hilbert spaces, spanned by basis vectors 
\x) and \x,/j) which satisfy the orthonormality conditions 



(x\x') — S(x, x') 



(x, fj,\x', u) = g liV (x)S(x, x'). 



(1.6) 
(1.7) 
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Strictly, we should here refer to GePfand triples [15] rather than abstract Hilbert spaces 
on their own, but we shall not be concerned with this standard of mathematical rigour. Of 
course, \x) is the familiar Dirac notation for the eigenfunctionals of the position operator 
which has continuous spectrum. They have a distributional nature and lead to a resolution 
of the identity in the form f d 3 x\x)(x\ = I. Moreover, the index of both \x,/j) and the 
associated 'bra' (x,/j,\ is viewed as that of a covariant vector density of weight \. The 
'Hamiltonian' operators H and H(a) associated to P and P Mi/ , respectively, are defined 
by (see (1.5) and (1.4)) 

(x\H \x) = P (x\x') (1.8) 

(x,/j\H(a)\x',v) = P*(a)(x,\\x',v). (1.9) 

The proper-time transformation kernels (also called 'heat kernels' by DeWitt [2,10], 
although in the mathematical literature one speaks about heat kernels for elliptic operators 
on manifolds endowed with a positive-definite metric) are defined by 

K${x^t) = (x,p\e- irH W\x',v) = K${T) (1.10) 

K (x, x'\ r) = {x\e~ iTHo \x') (1.11) 

where r is the proper-time parameter (not to be confused with the Euclidean-time param- 
eter of Euclidean field theory). The kernel K^,{r) (where the index n 'lives' at x and 
the index v 'lives' at x') is a solution of the initial-value problem consisting of the partial 
differential equation 

i^W = ^(«)^J(r) (1.12) 
subject to the initial condition 

K$(T = 0)= giu ,(x)6(x,xr). (1.13) 

For arbitrary values of a, the operator P^ v in (1.4), as well as P A = P px g Pil in (1.12), is 
non-minimal, in that the wavelike-operator part — g pv Q -\-R llu is spoiled by (l — ^) V u . 
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Nevertheless, if one knows K^l at a = 1, one can use this kernel, here denoted by K^J, (r), 
to evaluate K^),{t) according to the Endo formula [8] 




Equation (1.14) plays a key role in evaluating the regularized photon Green function (if 
no mass term or rotation of contour is used, one needs a regularization for Green functions 
as well, as we shall see following Endo [8]). Sections 2 and 3 perform this general analysis, 
while the full asymptotic expansion of the Feynman photon Green function at small space- 
time separation of the points x, x' is obtained in section 4. Coincidence limits of second 
derivatives of the associated Hadamard function are then evaluated in section 5, in light 
of their link with the regularized energy-momentum tensor. Concluding remarks and open 
problems are presented in section 6, while relevant details are given in the appendices. 

2. Photon Green functions in curved space-time 

As is well known, the photon Green function Gffi(x,x') = G °j, (unlike the contravariant 
realization in (1.1), we study hereafter the covariant one) satisfies the differential equation 



The formal solution of Eq. (2.1) admits, in the absence of zero and negatives modes, the 
integral representation (here g' = — det g^ v {x')) 



which is why we discussed the 'heat-kernel' K^J,(t) in the introduction. Bearing in mind 
that integration in (2.2) is taken over the positive half-line, we can say that it describes the 
massless limit of the Feynman propagator (for which one would have to add an infinitesimal 
negative imaginary mass). At this stage, the formula (2.2) needs a suitable regularization, 
because we use a Lorentzian-signature metric, so that contour rotation is not exploited 



^Wl^ =9^(x)S(x,x'). 



(2.1) 




(2.2) 
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to obtain a convergent integral. Following Endo, we use ^-function regularization and 
introduce a regularization parameter \ia defining (any suffix to denote regularization of 
the photon Green function is omitted for simplicity of notation) 

where Eq. (1.14) should be used on the right-hand side of (2.3). It should be stressed 
that the limit as s — > should be taken at the very end of all calculations, and cannot be 
brought within the integral (2.3) (see section 4). 

The kernel K^J,(r) is known as r — > and as a(x,x') — > through its Fock- 
Schwinger-DeWitt asymptotic expansion [2,13,14,16] 

oo 

where a = a(x, x') is half the square of the geodesic distance between x and x' , the bi-scalar 
\fK(x,x') is defined by the equation [10] 

A y/tf = <T.p V . (2.5) 

and the coefficient bi-vectors b n Mi y are evaluated by solving a recursion formula obtained 
upon insertion of (2.4) into Eq. (1.12). Such a recursion formula reads, for all n = 
0, 1, 2, oo, as 

v ,X K nu'-x + nb n ^ = (y/A 6n-i, M ^') A - R*b n - 1; xv'- (2.6) 

For example, one therefore finds bo ^> = g^ u ' , which is the parallel displacement matrix 
[2,17] along the geodesic between x and x'. The asymptotic expansion (2.4), called the 

local asymptotics of K^J,(t), only holds for r — > at small values of a(x,x'), which 
is precisely the framework of interest in renormalization theory and for applications to 
'laboratory physics'. At first sight, the above properties suggest therefore splitting the 
integral / Mi ,'(s, a) in (2.3) into an integral from to b plus an integral from b to oo, where 
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b e]0, 1[ to account for the r — > limit. However, as far as the asymptotic expansion of 
I^ v '{s^a) is concerned, we can use (2.3), (1.14) and integrate the local asymptotics (2.4) 
over the whole positive half-line of r, when x is very close to x' (see appendix A). Hence 
we find 

/•OO 

V(a,aO=/ r'^(r)dr~J^( a ) + J^( a ,a) (2.7) 



where 



1 oo 
~ T^M9' 1 J&nK i ' 1 " lf S ,n(^') (2-8) 



167T 2 

n=0 



~ 77^ X>V„V A (MAw^nliy^)) • (2-9) 

n=0 

In the asymptotic expansions (2.8) and (2.9) we have defined 

poo 

F s , n (x, x') = (a(x, x , )/2)'+»- 1 / y—^dy (2.10) 

./o 

and 

F fl>n (x, x'; a) = (a(x, x / )/2) T1 " 1 dr r s , j/" V»dy. (2.11) 

Jo J°°^ > 

At this stage, the Feynman photon Green function has the asymptotic expansion 

G$ ~<T* lim^^^W+J^^a))^)-! (2.12) 

where the asymptotic expansions (2.8) and (2.9) should be used, with integrals defined as 
in (2.10) and (2.11). 

3. Regularized integrals 

We are now going to evaluate the regularized integrals occurring in the space-time covariant 
form of the Feynman photon Green function (for more general results, see Refs. [18,19]). 
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For this purpose, we point out that the integral in the expression of F a ^ n (x, x') (see (2.10)) 
is a particular case of the integral 



/•OO 

I(/3) = / y-^e iy dy = iT(l - P)e~^. 
Jo 



(3.1) 



Recall now that the T-function T(z) = J °° y z ~ 1 e~ y dy, originally defined on the half-plane 
Re(z) > 0, can be analytically extended to a meromorphic function, with first-order poles 
at 0, —1, —2, — oo. With this understanding, we write that 

F a , n (x, x') = i(a(x, x , )/2) a+n - 1 T(l - 8 - n)e-^( s +") (3.2) 

where T(l — s — n) has first-order poles at 1 — s — n = —k, with k = 0, 1, 2, oo. 
To evaluate the double integral occurring in (2.11), we first exploit the identity 



a (x ,x' ) 
27 



/ " ,y- n e iy d y = i 3n+1 \r( 

I cz<t(x,x' ) \ 



r i -n, -i 



. cecr(x, x') 



2t 



-r 1 - n, -i 



. <r(x, x') 



2t 



where in square brackets we have the incomplete T-function 



(3.3) 



T(a, x) 



u a - x e- u du. 



(3.4) 



Hence we re-express F Sjn (x, x'\ a) in the form 



F s , n (x,x'-a) = (a(x,x')/2r + h^ +1 [l^'(a)-I^'(l)] (3.5) 



where 



TX, 

s, 



At this stage, we are led to consider the integral 



2r 



dr. 



/•OO 

J(J3,v,c)= / x^ _1 r(i/,ca;)da;. 



(3.6) 



(3.7) 



On setting y = cx, if Re(c) > 0, and exploiting the Leibniz rule and the fundamental 
theorem of calculus one finds [20] 



j= 1 



fid- 



because, for Re(/3) > and Re(/3 + z/) > 0, the total derivative of y^ u u ~ l e~ u du yields 
vanishing contribution to (3.8). 

We can however consider the analytic extension of T(/3 + v), after changing variable 
in the integral (3.6) according to ^ = T, which yields 

ee lim^T-^r (l - n, (e - i?^T) t) <£T (3.9) 

where a small positive £ > has been considered so as to be able to apply the result (3.8). 
In our case, (3 = —(s + 1), v = 1 — n, c = e — \ aa ^ x ) ; anc [ after making the analytic 
extension of Y{j3 + we find 

F s , n (x,*» = ^)/2) 5 +-i 3 ( s+ ") r( ( ~^~ ) n) (* s+1 " 1) (3-10) 

where r(— s — n) has first-order poles at s + n = k for all k = 0, 1, 2, oo. 

4. Asymptotic expansion of the Feynman photon Green function 

Our formulae (3.2) and (3.10) should be inserted into (2.8), (2.9) and (2.12) to work out the 
full asymptotic expansion of the Feynman photon Green function G^J,. For this purpose, 



it is crucial to take the limit as s — > in (2.12) at the last stage. Hence we find, as x 
approaches x' (which implies a(x,x') — > 0), 

where, after having defined 

u - a) s + ( i:'„ + ) 1 ( ;+ ) 1 ) v - vA < 4 ' 2 > 
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B nXul {s) = b n Xu/ y/A(x,x')(a(x,x')/2) s + n 



(4.3) 



we write 



g$ ( s ) = E T ^ - s - n ^ *) B n xAs)- 



(4.4) 



n=0 



What is crucial for us is the s — > limit of the sum (4.4). Indeed, on studying first, for 
simplicity, the case when the gauge-field operator reduces to a minimal (wavelike) operator 
(i.e. at a = 1), one finds 



(1) VAM^ 



n=0 



(4.5) 



having defined 



l\ /o\s+n 



fn nv' (s) = T(l - s - n)b n ^ (a(x, x )/2) 



Since 6q m^' = 5W we therefore find 



(4.6) 



(1) 2y/E( X ,x') VAM ^ 

v ' v ' n=l 



(4.7) 



which is very encouraging, since the first term on the right-hand side of (4.7) is precisely 
the first term in the Hadamard asymptotic expansion at small a(x,x') [16]. On the other 
hand, the Hadamard Green function is precisely the imaginary part of the Feynman Green 
function, in agreement with our formula (4.1). Eventually, we find therefore, at small 

a(x,x'), 



i \fK{x,x') 



8tt 2 a(x, x') + ie 



+ 



lim 



167T 2 s^O 

oo 

+ r(l - s - n)U n ^(s; a)B n Xv ,{s) 



n=l 



(4.8) 
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i.e. the 'flat' Feynman propagator, with +i£ term restored, plus corrections resulting from 
the gauge parameter (q / 1 leading to a non-minimal operator) and from non-vanishing 
curvature. 

A further crucial check is whether our infinite sum (4.4) is also able to recover the 
familiar \oga(x, x') singularity, which occurs for massive theories in flat space-time, and, 
more generally, even for massless theories (as is our case) but in curved space-time. For this 
purpose, it is enough to set a = 1 and focus on the sum in (4.5), having defined f n ^'(s) 
as in (4.6). Such a sum can be studied with the help of the Euler-Maclaurin formula (see 
appendix A and Ref. [21]), which provides, among the others, a term given by the integral 
(hereafter, since the discrete summation index n is replaced by the continuous variable z, 
we consider the coefficients b z ^ , functions of z that reduce to the coefficient bi- vectors 
b n ^ for z = n) 

vS^h = I °° r(1 - s - z)hz ^ (a(x ' x ')/ 2 ) 8+z ~ ldz 

= {a{x, x')/2) s T(l - s - z)b z ^e^- 1 ) ^ s (a(x,x')/2) dz 

+ ^ r(l - S - Z)b z ^e^" 1 ) ^{a{x,x')/2) dz ^ _ 



(4.9) 



At this stage we set z — 1 = y in the second integral in curly brackets on the right-hand 
side of (4.9), which therefore becomes 

V00 = (j^ + j^j r (~y ~ s)b y+1 ^ey lo ^^'Wdy. (4.10) 

We then begin to understand what happens: at small a(x,x'), the integrand in (4.10) 
becomes exponentially damped, so that the resulting asymptotic expansion of (4.10) is 
obtained from integration in the interval [0, y*] for some y* in a small neighbourhood of 
the origin. Here we first expand e y lo &{ a ( x ^ x ')/ 2 ) a t small y for fixed a{x, x'), and eventually 
take the a(x, x') — > limit. Such a procedure yields the non-uniform asymptotic expansion 

~ [V* 

Jpv'{s) ~ log(o-(x,a; / )/2) / yT(-y - s)b y+1 ^dy. (4.11) 

11 



On taking the s — > limit we therefore recover the familiar log(a(x, x')) singularity of the 
photon Green function, which results from non-vanishing Riemann curvature (in Minkowski 
space-time, the corresponding b y +i^ u > would instead vanish). 

5. Second derivatives of the Hadamard function in the coincidence limit 

In physical applications, one is interested in the energy-momentum tensor, which is ob- 
tained from the action functional (1.2) as 



rpflis _ 



2 SS 



Maxwell ' gauge ' 1 ghost 



(5.1) 



where (here a is the gauge parameter in (1.2), cf. Ref. [16]) 



rpfJ,U _ t?~1V _ 1 r 7 /3 T? 

1 Maxwell — F 7 4 r lf39 



(5.2) 



atT^ v 

gauge 



-A 13 A v — v A^ + 



A^ + \{A^ ^ 
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hi/ 



(5.3) 



rpflis 

ghost 



(5.4) 



On considering the Hadamard Green function, which here equals the imaginary part of 
the Feynman Green function, and is defined by 



> 



J + 



(5.5) 



one therefore finds, for the regularized energy-momentum tensor, the decomposition [16] 



fr£p,v\^ _ (7^) Maxwell + (T^) gauge + (T^ U ) ghost 



(5.6) 



where 



(T Ml/ ) Maxwell — T 1™ 
4 x'—>x 



x ^G^,. pT , + Gpy. Tp , — G^ T ,. p p, 

n H n H n H , riH , riH \ 

^rj'rfp' ^pP'-jr' <J /3p';r7' "t ^pr'wP' ^rp'-Pj'J 



(5.7) 
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a (TV) 



gauge 



lim 

x'^x 



-\g^E^^(Gl, lp ,-G%, pi ) 



(THghost = lim 



(5.8) 



(5.9) 



and we are assuming that our limits (5.7)-(5.9) do not depend on the choice of vacuum 
state. In the formula (5.8) we are using the DeWitt supermetric 



E ^u pr = g pp g ur + g pr g up _ g pu g pr 



and in (5.9) we consider the ghost Hadamard function [16] 



(5.10) 



G H (x,x') = ({x(x)^(x')] + ). 



(5.11) 



We should now specify in which order the various operations we rely upon are performed. 
Indeed, in the evaluation of the Feynman Green function in section 4, we first sum over n 
and then take the s — > limit. Here, we eventually obtain the energy-momentum tensor 
of the quantum theory according to the point-splitting procedure (5.6)-(5.9), with the 
understanding that the coincidence limit lim x /^ x is the last operation to be performed. In 
general, the analytic continuation in s and the coincidence limit do not commute [22], but 
the point-splitting result for can be made to agree with the local ^-function method, 
as has been proved in detail in Ref. [22] for scalar fields in curved space-time. 

It is clear from (5.7) and (5.8) that our analysis of the energy-momentum tensor 
is virtually completed if we can provide a closed expression for the coincidence limit 
linv^ x G^p,. ,. For this purpose we point out that (4.1)-(4.4) and the coincidence limits 
of appendix B show that the minimal-operator part of the Hadamard function contributes 
the divergent part (we write lim e ^ T( £ — k), with k = 0, 1, 2, in the formulae for such 

divergent contributions, where T(e — k) = j-^-j^r — I" 0(1)) 



lim 
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while the non-minimal operator part of the Hadamard function contributes further diver- 
gent terms given by (see (4.2)-(4.4)) 



(a — 1) lim 

e— >0 



r(^)(^/3 7 pT + Bp lpT ) + -T(e - l)C{3 lpT 



where (see appendix B) 

+ 

+ I 9\f3';pr' 





1 




9\f3'; 7pr' 


" 6 





R 



pr 



9xp'- X T ' 



Ryp 



9\f3>^p 



Rsyj- ~~\~ 



9\/3';yr' 



R p — 



9\f3';jp 



7pr' 



9\p 



(5.12) 



B 



(3-/PT 



4( 



°1 A/3'; t> 



9 ip' 



h x 



\9pr 



h X 

°lA/3'; 7 



+ 



a/3' 



9lT^j 

R\ 



12 



bl A/3' (R X T 9jp + -R A p fi , 7r)+ &1 r/3 



R~tp~\~ 



bi pp> R 1T ) 



+ 2 ( &1 P/3 ' ;7T ' ~ &1 T/9 ' ;7 '° ) 



6l A/3' (# A p 7T + ^ A r7p) + 75 



&1 



7/3' 



1? 



pr 



(5.13) 



pip- 



T = -2 ^ 62 T/3' 



#7P + 



>2 p/3' 



^77 



?2 7/3' 



9pr 



(5.14) 



We have therefore provided a covariant isolation of divergences resulting from every coinci- 
dence limit of second covariant derivatives of the Hadamard Green function. In manifolds 
without boundary, the work in Ref. [8] suggested that the trace anomaly resulting from 
the regularized T^ v has the coefficient of the Qi? term which depends on the gauge pa- 
rameter a (see, however, comments in section 6, end of second paragraph therein). In 
manifolds with boundary, the integration of such a total divergence does not vanish, and 
further boundary invariants contribute to the regularized energy-momentum tensor. Thus, 
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the calculation expressed by (5.12)-(5.14) is not of mere academic interest, but is going 
to prove especially useful when boundary effects are included (cf. results in Ref. [23]). Of 
course, physical predictions are expected to be independent of a, but the actual proof is 
then going to be hard. 

More precisely, the work by Brown and Ottewill [24] , which differs from our approach 
because the a = 1 case is there considered and the ^ and log (a) singularities in the 
propagator are there assumed rather than derived, has been exploited by Allen and Ot- 
tewill [25] to show that, on using the Ward identity and the ghost wave equation, the 
energy-momentum tensor is a-independent up to geometric terms, i.e. up to polynomial 
expressions of dimension length -4 formed from the metric, the Riemann tensor and its 
covariant derivatives. As far as we can see, our formulae (5.12)-(5.14) have precisely such 
a nature, having dimension* length -4 and being built from the metric and Riemann with 
its covariant derivatives (see appendix B), and have the merit of providing an explicit form 
of the general result of Allen and Ottewill [25]. Furthermore, the extension of these results 
to manifolds with boundary is, to our knowledge, an open research problem, and is of 
course relevant, for example, for the Casimir effect, which is a boundary effect in the first 
place. 

6. Concluding remarks 

By relying upon the regularized integrals (3.2) and (3.10), we have evaluated the asymptotic 
expansion of the Feynman photon Green function (see comments following (2.2)), here 
expressed in the form (4.1)-(4.4). Such an expansion corresponds to the singular part 
(i.e. divergent as a(x,x') — > 0) of the exact photon Green function [26,27]. We have 
endeavoured not to include mass terms because their addition 'by hand' spoils the gauge 
invariance of the original action nor is compatible with BRST invariance. Moreover, their 
addition to the Lagrangian of spinor electrodynamics leads to a photon propagator with 
a k° part in momentum space, incompatible with perturbative renormalizability; this is 

* the coincidence limits of b\ ^ v > and 62 \xv> have dimension length -2 and length -4 , 
respectively, each covariant derivative has dimension length -1 , while Riemann, Ricci and 
the scalar curvature each have dimension length -2 . 
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compensated by adding an auxiliary vector field which spoils unitarity [28] . It is therefore 
rather important to study photon Green functions from the point of view of massless 

2 

QED theory. The calculational techniques used in the presence of a ^A^A^ term in the 
Lagrangian [2,10,16] are then no longer available, which is why our formulae (4.2)-(4.4) 
provide a novel way of expressing the familiar singularities in the asymptotic expansion of 
the Feynman photon Green function. 

Our second original contribution is a covariant isolation of divergences in the energy- 
momentum tensor resulting from every coincidence limit of second covariant derivatives 
of the Hadamard Green function, including the contribution of the gauge parameter a 
which leads to a non-minimal operator on the potential in the path integral. In Minkowski 
space-time it is well-known how to deal with arbitrary a in the photon propagator [5], 
and hence it is desirable to deal with arbitrary a also in curved space-time. One then 
finds, as we have done in (5.12)-(5.14), an explicit expression of the geometric terms of 
dimension length -4 up to which the regularized energy-momentum tensor is known to 
be a-independent, according to the general analysis of Allen and Ottewill [25]. The a- 
dependence of the electromagnetic trace anomaly found in Ref. [8] is instead an incorrect 
claim, since the analysis of Nielsen and van Nieuwenhuizen [11] has shown that the corre- 
sponding a-dependence of the ghost effective action cancels such an a-dependence [11,25]. 

After the early work in Refs. [7,8,26,29], for example, there has been recent work by 
other authors. More precisely, the work in Ref. [18] has obtained the Euclidean Green 
function for an operator of Laplace type (corresponding to the choice a = 1), while the work 
in Ref. [19] has obtained results of very general nature, including in particular a recursive 
algorithm which holds for non-minimal operators (a ^ 1 in our case) with positive-definite 
metrics. Heat-kernel asymptotics for non-minimal operators has instead been worked out 
in detail by Gusynin and his collaborators (see, for example, the work in Refs. [30,31]). 

We should admit that we might have based all our analysis of Green functions on Eq. 
(2.31) of Ref. [32], but our approach in sections 2-4 is more directly related to the properties 
resulting from the Fock-Schwinger-DeWitt asymptotic expansion in the case of minimal 
operators. Casimir energies will eventually motivate the use of a modified Schwinger- 
DeWitt ansatz for manifolds with boundary [33], along the lines of the work in Ref. [34]. 



16 



The central object in Ref. [34] is indeed the scalar Feynman Green function, which is a 
first step towards the photon Green function in curved space-time, eventually including 
boundary surfaces, which are crucial in the Casimir effect [35,36]. At that stage, our 
formulae (5.12)-(5.14) will be part of the general formula leading to the regularized energy- 
momentum tensor. Hopefully, such investigations will find application to the exciting new 
problem of Casimir apparatuses in a (weak) gravitational field [34,37], where the regularized 
energy-momentum tensor of QED in curved space-time can be used to evaluate the force 
acting on the Casimir apparatus. 
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Appendix A 

We are now aiming to discuss under which conditions one can use (2.3), (1.14) and integrate 
the local asymptotics (2.4) over the whole positive half-line of r, relying upon private 
correspondence with Ivan Avramidi. For this purpose, we start with the better defined 
Euclidean theory, where the heat kernel [38] has truly such a nature, and has the following 
asymptotic expansion in four dimensions: 



where Pn, of a polynomial nature, denotes the first N terms and Rn is the exponentially 
small remainder, for which 



U(x,x';t) = (4nt)~ 2 e~ 



a (x , x 1 ) 
2t 



[P N (x,x';t) + R N (x,x';t)} 



(Al) 




(A2) 



Now by temporary addition of a finite mass term, the integral 




(A3) 
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becomes meaningful (strictly, for m such that the second-order differential operator ruling 
the field becomes positive-definite), and the contributions of Pn and Rn can be integrated 
separately. At this stage, the mathematical limit as m — > 0, and Wick rotation back to 
'real time', make it possible to obtain the asymptotic expansion (2.10). 

The Euler-Maclaurin formula [21] used in section 4 asserts that, if / : [0, oo[^ R is a 
function having even-order derivatives which are absolutely integrable on (0, oo), then, for 
all k = 1, 2, oo, 



where Bi s are the Bernoulli numbers, and R m (k) is the remainder term, majorized by 



Appendix B 

In the course of deriving and further elaborating the formulae (5.12)-(5.14), we need the 
following coincidence limits [39,40], here denoted by square brackets [...] as in 





m — 1 






Refs. [2,10,17]: 



M = W = M = o 



(Bl) 




(B2) 




] — W;Xnu'] — k;Xjuv] — k;A>V 




(S3) 




(B4) 



(B5) 




= 5 



(B6) 
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//d'':nj — ypu';a'\ — 



[g^u'-ap] — -[g/j,v';af3'] — --^Riiuap 

[gpu';ctP-y] = ~ g (yRjjivafi;^ + Rp.uct-f.p^j 

{gnv';af3~f5] + [g p.v' -pct-ys] + [g^u'-'yaps] + [g nv' ;5ct/3j\ 

+ g (R X /3aj + R^afij Rpv\5 + {R X p a S + ^'ja/?) Rpu\-y 

[fi , M^';a/37<5] ~~ [g^u'-.^a^s] = — Rpuap-fS ~ -^R^sR^XaP 
[g^u'-oip-fS'] = —[g^v'-aP-yS] + \g f iu';al3'y];5 



[VA] = 1 [a/A J = [VA ;M ,] = 

[A/A;,,] = \R,U = "[a/A;^]- 



;a/37 
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{RaP;y + Ra-/;P + -R/3 7 ;a) 



otp^jb 



pctp'yS 



- g [R ap R p p lS + Rp P Rp c 

+ R 1P R p a ps + R8pR p a [3^ + g { l Ra P S p f3l5 
+ RppS p aiS + Rj P S p a/3S 

+ Rs P s p aPl ) 



-(R a pR~/5 + Ra-yRps + RasRpj) 



aPjS' 



ctp^/S 



+ 



a/37 



■5 



bi pv> 



Rg^v Rpv 
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>2 fiu' 



1 1 1 1 

— qRR(j,v + g I I Rgv + 7^RppR P v ~ P ii^^o pv 



+ {j2 R2 + 3^ QjR " 1^^ + 1^^^) (S19) 



&i ^';p 



R- p + 2R \\jg^i, R^u; P ^R^p-^p 



(520) 



b\pLv';pLU 



1 



1 



i r i 

+ g | ~~ 7^R9tivRpuj — R\xv, 



puj 



+ --Rf, Rxvpcu — —RR 



l 



l 



12 



pVpUJ 



R pRp, v \uj R ^RpvXp 



+ 9^9^ 



6 



9 nv' pu 



6 



(521) 



+ 



(522) 



As is stressed in Ref. [17], the possibility of taking such coincidence limits relies on the 
assumption that the space-time points x and x' have a unique geodesic passing through 
them. This is the case if the points are close enough to one another, but there are physical 
instances where it does not hold. The world function is then no longer single-valued, and 
the existence of partial (or covariant) derivatives is not obvious a priori. A global theory 
of the world function a, covering such singular cases, is very complicated and goes beyond 
our aims; we therefore always assume that the geodesic passing through x and x' is unique 
and that partial derivatives exist [17]. 
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